An integral equation for distorted wave amplitudes 
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We derive a new integral equation that allows the calculation of the scattering or annihilation 
amplitude of two particles subjected to two potentials when the corresponding amplitude for one 
potential only is known. We assume that scattering or annihilation occurs through one of the 
potentials, while the other potential affects the particle wave functions. Our expression is valid for 
any choice of the distorting potential and for any particle model. Our technique does not require 
the expansion of the amplitude into partial waves, and allows the study of models that are generally 
difficult to solve by means of the Schroedinger equation. 
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There are many examples in nature in which a physical 
system is affected by two distinct potentials. For exam- 
ple, nuclear beta decay occurs through weak interactions, 
but electromagnetic forces between the charged nucleus 
and the emitted electron affect the motion of the final 
particles by decelerating them [2 |2] . Another example 
is nucleon-nucleon scattering [3], in which the exchange 
of virtual mesons before the scattering can be described 
by an effective potential that enhances the amplitude for 
the process. A similar situation has been invoked in dark 
matter studies to generate anomalously high annihilation 
cross sections at low velocities in relation to recent cos- 
mic ray data (an effect dubbed the Sommerfeld enhance- 
ment after Ref. [4j). Another example is provided by the 
quark-antiquark bound state system [5] , in which there is 
a potential for quark confinement and another for quark 
decay. Still another example occurs in the minimal super- 
symmetric standard model when the lightest and second 
lightest supersymmetric particles are close in mass and 
convert into each other prior to their final annihilation 

m- 

In all of these examples, one potential (from here on 
indicated with W) is responsible for the final scattering, 
annihilation or decay, while another potential (here re- 
ferred to as the distorting potential V) affects the initial 
particle wave functions and thus modifies the transition 
amplitude. 

In this Letter, we derive an exact integral equation 
for the transition amplitude in the presence of two po- 
tentials V and W in terms of the transition amplitude 
with only one of the potentials, say W. We also de- 
rive the equation relating the corresponding initial-state 
wave functions for two-particle states. To the extent of 
our knowledge, these equations have not been presented 
before in the literature. 

In our derivation, we do not make any assumption re- 
garding the strength of one potential compared to the 



other. In other words, we do not invoke the so-called 
Distorted Wave Born Approximation (DWBA) [3 H] , al- 
though we discuss the simplification of our integral equa- 
tion when the DWBA is applied. 

Although our derivation is general, we present it here 
for the common case of two incident particles. The fi- 
nal state is completely arbitrary. Consider a process in 
which two particles xi s-nd X2, of masses mi and 7712, re- 
spectively, scatter or annihilate into a multi-particle final 
state. 
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The particles that take part in the process can be 
fermions or bosons. In the center-of-mass (CM) frame, xi 
has momentum pi = p, and X2 has momentum p2 = —p. 
If one or both particles have spin, we quantize the total 
spin of the initial two-particle state along the axis defined 
by the direction of pi. Thus, if Ai and A2 denote the he- 
licities of particles 1 and 2, i.e. their spin projection in 
the direction of their respective momenta, the total CM 
spin of the initial state has projection Ap = Ai — A2 along 
the direction of p. 

In labeling the initial two-particle state we need to dis- 
tinguish between plane-wave states, scattering waves in 
the presence of the potential W alone, and scattering 
waves in the presence of both potentials V and W. We 
denote these states by \pSApT]), \ipps%r,)' and ItppSA^r)^ 
respectively. Here S is the total CM spin, the index ± 
distinguishes between incoming (— ) and outgoing (-1-) 
scattering waves, and 77 (later suppressed) denotes ad- 
ditional quantum numbers besides momentum and spin. 
We also use a condensed notation a = {p,S,Ap,i]) for 
the ensemble of quantum numbers. 

We remark at this point that for the process in Eq. (fTl) 
with xi and X2 in the initial state, one should use out- 
going scattering waves since it is these waves that reduce 
to plane waves in the infinite past {t — >■ — c»). For this 
reason, in the following we restrict ourselves to outgoing 
(-I-) waves. 

The states IV'^^) satisfy the Schroedinger equation 



iHo + W)\yj^+)=E^\^^+) 



(2) 



Here, Hq is the free hamiltonian, and according to the 
prescription in scattering theory, Ea is the total energy 
of the two particles in the CM frame, equal to the total 
energy in the absence of the potential: 



HqIu) = Ea\a). 



(3) 



In the non-rclativistic limit, Ea = p^/(2/i), where ^ — 
mim2/{mi + TO2) is the reduced mass of the initial par- 
ticles. 

The states \'ip'^~^) satisfy the completeness relation 
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Eq. (l7| becomes 



\^i) = \<^)+ Id, ^'^^^^'''"' 



Ea — Ey + ie 



(10) 



The quantity T^q, which coincides with the difference 
of the T-matrices for the process in Eq. (IT]) when both 
V and W or only W are present, satisfies the Lippmann- 
Schwinger equation 



T - V^ 



di 



Ty^, (11) 



Ea — Eyi + ie 



where 



for an appropriate integration measure that can contain 
a sum over discrete quantum numbers. The states \tp^~^) 
are normalized so that 



{r^^ + \i^y)^{2nrS{a-a'), 



(5) 



where similarly to the integration measure, the Dirac 
delta function can include Kronecker deltas in the case of 
discrete eigenvalues. We remark that the symbols a and 
a' appearing in Eqs. Q and ([5]) contain both continuum 
and bound states. 

We now add the distorting potential V. The states 
\ijj'^) satisfy the equation 



{Ho + W + V)\r^+)^E^\^+), 



(6) 



where, according to scattering theory, E^ is the same en- 
ergy appearing in Eq. pi). Following standard procedure 
(see, e.g., [ZIIH]), Eq. ([6f^can be cast into the Lippmann- 
Schwinger form 

\^i) = \^^+) + {E^-H,-W + ter'V\i,+). (7) 

To see that Eq. ^ and Eq. ^ are identical, one can 
multiply Eq. (M) on the left by the operator Ea — Hq — 
W, after which the term {Ea — Hq — W)\%1^^^) vanishes 
because of Eq. ([2|. The extra term +ie in Eq. ([7| is 
introduced in order to deal with the singularity in the 
operator {Ea — Hq — W)~^ , and takes into account the 
sign (-I-) describing the out-going state \ipa)- 

From Eq. ^ we derive our main result as follows. We 
insert a completeness relation^ffl) into the second term 
on the right hand side of Eq. (M), thus 

{Ea-HQ-W + l€y^V\ll^t) = 

d,{Ea~HQ-W + ie)-'\^'^+){i^'^+\V\i^t) = 



d'y 
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Ea ~ E~, + ie 
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The label 7 describes the set of quantum numbers asso- 
ciated with the intermediate state \tjj^^)- Defining 



T,a = {i^^+\V\iji), 



(9) 



yj^, ^ (vr+|F|^!r+) 



(12) 



are the matrix elements of V in the basis \ip^^). Eq. ( 11 ) 
can be obtained by multiplying Eq. ([7]) on the left by 
('(/'^''^l V, and inserting a completeness relation in a way 



similar to the derivation of Eq. ( 10 ) 



We now use Eg . (11) to eliminate T^a from Eq. (10 1 
We expand Eq. ( 11 ) recursively in powers of V^ , 
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and insert this expansion into Eq. (10 1, 
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{Ea -E^+ ie){Ea ~ Ey + ie) 



(14) 



We finally notice that the last expression corresponds to 
the series expansion of 



IV'^) = ie+) + / d7 ''^^^^^'^" 



Ea ~ Ej + ie 



(15) 



Eq. (151 is one of our two main results (the other one is 
Eq. (18 I below). 

Eq. (15) relates the state |V'a), solution of the 



Schroedinger equation in the presence of both poten- 
tials V and W, to the state IV'^'*'); solution of the 
Schroedinger equ ation in the presence of the potential 
W only. Eq. (15) can be used to find the relative nor- 
malization of the Itfja) and IV'^'*') states. 

We stress that no assumption regarding the strength of 
V relative to W has been made in deriving Eq. ( 15 ). The 



potential W does not enter Eq. (15) directly, so to find 
IV'q) one only needs to know the states \ip^'^) and the 
potential V. We remark that Eq. ( 15 ) cannot be obtained 
from Eq. ([7| by simply inserting a completeness relation 
in the way we derived Eqs. (10) and (11). 



From Eq. ( |15[ ) we derive an integral equation for the 
transition amplitude of process M. Introducing the fi- 
nal state 10^), where /? labels the quantum numbers of 
the final particles x'li ■••jXn; the transition amplitude for 
process il\ when only W is present is 



Ml^ = (0^|V^ 



H'^ 



(16) 



We refer to A^a^, as the undistorted amplitude, undis- 

torted by the potential V. When both V and W are 
present, the transition amplitude is 
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We refer to M'l^ as the distorted amplitude. Multiplying 



Eq. (15) on the left by (^^g], we obtain 



■^^a 



■M^a 
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Mt V^ 

P^ lot 

E„- E^+7 



(18) 



This is our second main result. Eq. (18) is an integral 
equation for the distorted amplitude, which can be solved 
if the undistorted amplitude and V^ are known. 



In Ref. [TT] we show that Eq. ( 18 1 is the non-relativistic 
limit of sums of ladder diagrams in quantum field theory. 
Here we only sketch the proof graphically in Fig. 1. In 
this case, V and W are effective potentials that may in- 
clude radiative corrections to any number of loops. 

The graphical derivation allows the understanding of 
which quantum field theory diagrams are included in our 
Eq. (18). They are ladder diagrams containing multi- 
ple powers of the effective potential V, while the other 
effective potential W, responsible for the final inter- 
action, appears only once. Radiative corrections like 
bremsstrahlung of particles from the initial or intermedi- 
ate legs can also be treated in our equation by including 
the bremsstrahlung particles in a modified final state. 
Notice in this regard that non-relativistic quantum me- 
chanics does not allow for the number of particles in the 
final state to change, so an ni-particle process with extra 
712 bremsstrahlung particles should be considered as an 
ni-|-n2-particle process. 

So far our derivation of Eqs. (15) and (18) has been 
exact. We now discuss Eq. (18) in the Born approxima- 



tion with respect to the potential W. In the first Born 
approximation, 



l< 



W+\ 



OiW), 



(19) 



where \a) is a free state of Eq. ( [3|). To the same order in 
W, the matrix elements in Eq. (^12| become 



V^Z^{l\V\a)+0{W). 



(20) 



Thus the distorted amplitude Mt^ in the first Born ap- 
proximation in W, which is known as the DWBA ampli- 
tude, satisfies the integral equation 
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(21) 










FIG. 1. Graphical derivation of the ladder equation that re- 




duces to our Eq. ( 18 1 in the non-relativistic limit. Gray circles 



stand for the effective potential W , including radiative cor- 
rections to any number of loops. White ovals indicate the 
effective potential V to any nuniber of loops. The diagram on 
the LHS, indicating Mt^ [Eq. (|l7|], is the sum of an infinite 

ladder of effective potentials V and a final W interaction. The 
first diagram on the RHS stands for A^Jj^ [Eq. (16l]. 



The only difference from Eq. (18), besides the Born ap- 
proximation implicit in A^a^, is the presence of V^a = 

(7|V^|a) instead of V^^^^ — {i^^'^\V\il)^,'^), i.e. plane wave 

states replace W scattering states in the matrix elements 
of V. Eq. (21) is a good approximation whenever the 



potential W is weak, which is a common case. 

For potentials V that do not depend on the spin of 
the particles, Eq. (21) with the Born approximation in 



W splits into separate independent equations for each 
value of the total spin S. (If W is also spin-independent, 
then an analogous derivation shows that the full Eq. ( 18 ) 



separates into a set of independent equations for each 
total spin 5*.) To see this, we factorize the two-particle 
free state \a) = jpS'Ap?]) into an orbital part jp?/) and a 
spin part jS'Ap). Thus, 

\a) ^ jp^ApTy) = |p77) \SK^). (22) 

Similarly, we factorize the intermediate state into 

|7)^|k^X'7') = IW)|5'A;,). (23) 

Here k is the internal three-momentum. Notice that for 
5 = 0, only the orbital parts exist: \a) — jpyy) and [7) = 

IW). 

For simplicity, we define the -f z axis along the initial 
momentum, and specify the direction of the internal mo- 
mentum k using spherical coordinates (0^, (j)}S). Then the 
helicity state |S"AJ^) with spin polarization along k can 
be obtained through a rotation from the helicity state 
|S"A' ) polarized along p. In the convention of Jacob and 
Wick (see Eq. (6) in Ref. [3]), the two helicity states are 
related through a rotation of Euler angles (0fc, ^fc, — 0fe), 

|5A;,)=7^(0,,0fc,-(/.fc)|5A;,). (24) 

The rotation operator TZ{<j)k,&k,-'4>k) is related to the 

( S) 
Wigner Z?-matrix Dj^^^j, through 

(5Ap|5A;,> = (5Ap|7^(</>fc,0fc,-0fc)|5A' ) = 



D 



{S) 



-<pk) 



=,-'(Ap-A;)0fc AS) 



"■A^.A' 



(25) 



{ S) 

The function rf)vfM'(/^)' which is a real function and is 
known in the literature as the small d-matrix, is defined 
using the generator Jy of rotations about the y-axis be- 
tween spin states polarized along the z axis: 



d^-il^) = {SM\e-^'^^/'-^\SM') 



Using Eqs. p'A 



jugate of Eq. ( 



(23), and the complex con- 



25 ) with D 






c, -</>*:) 



-^A A' {~4'k,0k,4>k), the matrix elements V-^a — {l\V\c() 
of a spin-independent potential factorize as 

V,^ = (5XI(k|V^|p)|^Ap) = {S'A'^lSAp) (k|1/|p) 
= ^ss' D[^Jj,,^i-^k,Ok,^k) (k|l/|p). (27) 

Here (k|l/|p) is the Fourier transform of the distorting 
potential V{r), 

(k|T>|p) = f d^rVir)e-'^^-P^-'' = V{p-k). (28) 



Hence 



l(S) 



V^a. = Sss'D)^^Jj^,^i-cj)k.,ek, <Pk) ^(P - k). 



(29) 



Notice that for spin zero particles, one sim ply obtains 
V^a = V"(p-k). The Kronecker delta in Eqj29| forbids 



any mixing of different total spins in Eq. (21 ). 



Eq. (21 1 then splits into a separate set of 25 -f 1 equa- 



tions for every value of the total spin S. To display the 
structure of the each set clearly, we introduce the abbre- 
viated notation 

MUv)^Ml^SK^, M%dv)^M%sK,^. (30) 

Then for A = -S", -5 + 1, . . . , 5, we find 

>'sa(p)--^5a(p)+ 

^^J (2^)3 Ep-Ek+ie 

(31) 



In general, Eq. (31 1 has to be solved numerically for each 



value of S to obtain M.g^{p) 



Eq. (31 1, obtained in the DWBA, is valid for any spin 
of the annihilating particles. In particular, for bosonic 
(26) particles or for particles which are in a singlet state, S = 



A = 0,<],(- 



M+oiP) ^ M'^ooip) 



j) = 1, and 



d^k V{p- 



k)X+(k) 



(27r)3 Ep - Ek 



(32) 



This singlet equation, Eq. (32 1, appears in Ref. [10] in a 



form limited to the specific case of a Yukawa potential 
V{r). Our equations apply to a generic potential and to 
any value of the total spin. 

W e h ave derived a new exact integral equation, 
Eq. (18), that gives the amplitude Mt^ for two parti- 
cles to scatter or annihilate into a multi-particle state in 
the presence of two potentials V and W once the am- 
plitude A^ a„ for the same process in the absence of the 

potential V is known. We have also obtained a similar 
relation connecting the ini tial particle states in presence 
and absence of V , see Eq. ( 15 ). We have also shown how 
our equation simplifies in the DWBA, Eq.(21 1. For spin- 



independent potentials V , our DWBA equation separates 
into independent sets, Eq. (31 ), one for each value of the 
total spin S of the co llid ing particles. The same is true 
for our full equation ([18) if also W is spin-independent. 
Our integral equations ( |18[ ) or (31 1 can be solved numer- 
ically for the amplitudes Mt^ or 7V{J^(p), respectively, 
without an expansion in partial waves. Avoiding the par- 
tial wave expansion is convenient, for example, when the 
s-wave is slightly suppressed or there are p-wave reso- 
Numerical studies of Eq. (31 1 will be presented 

ID. 



nances 

in a separate paper 



This work was supported in part by NSF award PHY- 
0456825 and NASA contract NNX09AT70G. 



[1] E. Fermi, Nuovo Cim. 11, 1 (1934); Z. Phys. 88, 161 

(1934). 
[2] V. B. Berestetskii, E. M. Lifshitz and L. P. Pitaevskii, 

Relativistic Quantum Theory, Pergamon Press (1971). 
[3] S. M. Dancoff, Phys. Rev. 78, 382 (1950); E. E. Salpeter 

and H. A. Bethe, Phys. Rev. 84, 1232 (1951). 
[4] A. Sommerfeld, Ann. Phys. 11, 257 (1931). 
[5] R. Delbourgo, A.Salam and J. Strathdee, Nuovo Cim. 

ser. 10, 50, 193 (1967). 
[6] J. Hisano, S. Matsumoto, M. M. Nojiri, O. Saito, Phys. 



Rev. D 71, 063528 (2005). 
[7] A. S. Davydov, Quantum mechanics, Pergamon Press, 

London (1965). 
[8] S. Weinberg, The Quantum Theory of Fields - Founda- 
tions, Cambridge University Press (1995). 
[9] M. Jacob and G. C. Wick, Ann. Phys. 7, 404 (1959). 
[10] R. lengo, JHEP 0905, 024 (2009). 
[11] L. Visinelli and P. Gondolo, in preparation. 



